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We address the problem of a single "spin-down" impurity atom interacting attractively with a 
spin-up Fermi gas in two dimensions (2D). We consider the case where the mass of the impurity is 
greater than or equal to the mass of a spin-up fermion. Using a variational approach, we resolve 
the questions raised by previous studies and we show that there is, in fact, a transition between 
polaron and molecule (dimer) ground states in 2D. For the molecule state, we use a variational wave 
function with a single particle-hole excitation on the Fermi sea and we find that its energy matches 
that of the exact solution in the limit of infinite impurity mass. Thus, we expect the variational 
approach to provide a reliable tool for investigating 2D systems. 



The impressive realization of Fermi systems in ultra- 
cold atomic gases has greatly renewed interest in funda- 
mental problems in pairing phenomena. For attractive 
interactions, the binding of fermions can alter the statis- 
tics of low-energy excitations and thus determine the low- 
energy behavior of Fermi systems such as Bose-Fermi 
mixtures and two-component Fermi gases. This phe- 
nomenon in Fermi systems is nicely captured by the fun- 
damental problem: What is the ground state of a single 
impurity atom interacting attractively with a Fermi gas? 
This problem differs from the usual polaron problem of an 
impurity particle interacting with a background medium, 
since the medium is fermionic rather than bosonic. It 
also has connections with other Fermi systems: For a 
fermionic impurity, it corresponds to the extreme limit 
of spin imbalance in a two-component Fermi gas • 

In three dimensions (3D), it is known that this "spin- 
down" impurity can undergo a sharp transition to a new 
ground state by binding a fermion from the spin-up Fermi 
sea 0,0]. For weak interactions, the impurity is initially 
dressed with density fluctuations of the Fermi gas, form- 
ing the so-called "polaron" state 0, 0] ■ Then, with in- 
creasing attraction, the impurity eventually binds a spin- 
up fermion to form a dimer or molecule state in the case 
of equal masses (m^ = m\) 0-ll|. For a sufficiently light 
impurity (m^ < m^), there is even the prospect of it bind- 
ing two spin-up fermions to form a dressed trimer 12 1 . 

However, the situation is less clear in lower dimensions, 
where quantum fluctuations are important. For 2D Fermi 
gases, Ref. [l3| find no evidence of a polaron-molecule 
transition when particle-hole excitations of the Fermi sea 
are neglected in the molecule state. Instead, they find 
that the polaronic description persists into the regime 
of strong attraction, even when it is no longer accurate. 
This suggests that an extension of their variational ansatz 
is required [l3| . In this paper, we resolve this ques- 
tion and show that there is indeed a polaron-molecule 
transition in 2D, provided one treats the polaron and 
molecule on an equal footing and includes particle-hole 
excitations of the Fermi sea in both variational wave func- 
tions. Here we consider a molecule wave function with 
a single particle-hole pair and we show that this is suffi- 



cient to reproduce the exact impurity energy in the limit 
of infinite impurity mass. We also show that the vari- 
ational approach can never produce a polaron-molecule 
transition in ID, even with the inclusion of particle-hole 
excitations, which is consistent with the finding that the 
polaronic description appears accurate across the whole 
range of interactions in ID [l^. We expect our results to 
be important for the ongoing investigation of 2D atomic 
Fermi gases [3, [lB| , as well as for other 2D Fermi systems 
in condensed matter such as electron-hole bilayers 

Model In the following, we consider a two-component 
(t, i) 2D system with short-range interspecies interac- 
tions, described by the Hamiltonian: 
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where eko- = (we set h = 1), 57 is the system area 
and g is the attractive contact interaction. For a Fermi 
mixture, the Hamiltonian ([T|) completely describes the 
low-energy behavior, since Pauli exclusion suppresses s- 
wave interactions between the same species of fermion. 
For bosonic spin-down particles, we must also include in- 
traspecies interactions between bosons. However, Eq. ([T]) 
is sufficient for this work since we restrict ourselves to the 
problem of a single J,-particle in a f Fermi sea. 

To make our results independent of the UV cut-off for 
the short-range interaction, we use the energy eigenvalue 
equation for the two-body (t and J,) problem: 
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where A is the momentum cut-off and £b is the two-body 
binding energy, which always exists for an attractive in- 
teraction in 2D. This allows us to replace g with eb h^ 
our calculations and then take A — >■ co. 

We parameterize the non-interacting Fermi gas with 
the Fermi energy sp = fc|,/2m-|-, where kp corresponds 
to the Fermi wave vector of the spin-up Fermi sea. Defin- 
ing the mass ratio r ~ m^/m^, we focus on the regime 
r < 1; the limit ?' — > can be compared with the exact 
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FIG. 1: (Color online) Energy of an impurity with infinite 
mass in the strong-coupling limit eb/sf S> 1. The thick 
(red) line corresponds to the exact solution ((HJ, while the 
thin (black) line with data points is the energy of the dressed 
molecule calculated using the variational wave function in 
Eq. ((S]). The energy calculated using the polaron varia- 
tional wave function ([3]) is depicted in the inset. In the 
limit eb — s- oo, the polaron energy approaches E + eb ~ 
—Ep/ie — 1) ~ — 0.582eF. Note that E + eb depends linearly 
on 1/eb for the polaron in this limit, which gives rise to the 
apparent kink in the inset for the units 1/ \n{EB /ef)- 



solution [13[ , while the equal- mass case r = 1 corresponds 
to current Fermi gas experiments [l^. 

Polaron For weak interactions (or, equivalently, large 
densities) where Eb/ef 1 , the ground state is expected 
to be a polaron with the approximate wave function 
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which contains just one particle-hole excitation on top 
of the Fermi sea |A^) with N spin- up fermions. We 
determine the ground state energy E by minimizing 
(P| {H — E) \P) with respect to the amplitudes a, yield- 
ing the implicit equation: 
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where E'kq = —E + fq-ki + Ckf — fqt- The real eigen- 
state will, of course, include terms with greater num- 
bers of particle-hole pairs, but it has been argued that 
these terms approximately cancel in the variational equa- 
tions for E [8|. Indeed, if we take akq to be indepen- 
dent of q, i.e. akq = /3k, then we obtain a closed set 
of equations for ao, /3k that do not involve higher-order 
terms. Thus, we can estimate the "exactness" of our 
variational approach by comparing the energy of the full 
variational wave function with that of our q approxi- 
mation, where we simply replace i?kq in Eq. ^ with 
(i^kq)q = -E + ek^ + ekt + Spir - l)/2. Note that, for 



r = 1, this is equivalent to replacing ii'kq with i?kO: but 
it provides a better estimate of the energy for r < 1. 
We find that the q approximation becomes better with 
decreasing eb/ef, implying that |P) is accurate at weak- 
coupling, as expected. 

Dressed molecule For stronger interactions, we con- 
sider a molecule wave function with a single particle-hole 
excitation of the Fermi sea: 
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Note that we use a Fermi sea of — 1 spin-up fermions 
in order to preserve particle number. By minimizing 
(Af I {H — E) \ M) with respect to the (/s's, we obtain the 
simultaneous integral equations: 
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where Em = E + ef, i?kk'q = -Em + Cq-k-k'i + Ckf + 
Sk't - eqt and G(k,q) = ^ J^w 'Pkk'q/Ek' <^k'- Here, it 
is understood that the hole momenta q satisfy |q| < kp 
and particle momenta k satisfy |k| > kp. For G'(k, q) = 
0; Eq. ^ recovers the equation for the "bare" molecule 
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considered in [13| . 

with center-of-mass momentum Q = 0. The bare state 
|Mo) never has a lower energy than |P) and is unable 
to reproduce the strong-coupling limit for r = 0. In this 
case, the impurity energy is given exactly by [l3| : 



'Ze p 

E = —Eb / dkk5[){kkp) 



(8) 



where cot{So{kkp)) = \n{k'^ e p / e b) / ■ For eb — >■ oo, we 
have Jo — and thus E ~ —eb—Ef, which disagrees with 
the bare molecule energy oi —eb- However, the inclusion 
of G(k, q) has a dramatic effect on the molecule energy. 
To sec this, we first assume G(k, q) is small and drop 
the sums involving G'(k, q) in Eq. This allows us 

to solve for G(k, q) and then substitute it into Eq. 
For r < 1, the solution of the resulting implicit equation 
yields E + ef + Eb oc — In(eB), a correction which clearly 
diverges for eb — > oo. This implies that G'(k, q) is not 
perturbative in and thus cannot be neglected in the 
strong-coupling regime. 

We have solved Eqs. I©-© numerically by discrctiz- 
ing momentum space and converting the problem into a 
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matrix equation. The results for an impurity with infi- 
nite mass (r = 0) are plotted in Fig. [1] We see that the 
energy curve for \M) collapses onto the exact result as 
we take Eb ^ oo, thus demonstrating that a molecule 
wave function with a single particle-hole pair is suffi- 
cient for capturing the strong-coupling limit. By con- 
trast, the energy for the polaron \P) rapidly approaches 
E+eb = —EFl[e—l) when eb — > oo, as depicted in Fig.[T] 
(inset). This limit is obtained analytically by assuming 
E + Eb ^ Eb and approximating the expression in the 
q-sum in Eq. (g]) to get E ~ -eb ln(l - 2ef/{E + Eb))- 
Taking exponentials and expanding in {E + eb)/eb then 
gives the final result. For general r < 1, we can test the 
accuracy of |M) by comparing its energy with that of the 
q approximation, where we replace -Bkk'q in Eq. ([7]) with 
(-Bkk'q>q = -£'j\/+ek-i-k'4,+ekt+ek't+^F(f-l)/2. As an- 
ticipated, wc obtain excellent agreement for eb/ef 3> 1. 

Polaron-molecule transition The fact that the po- 
laron state \P) has a much higher energy than the 
molecule state |M) in Fig. [1] already suggests that there 
must be a polaron-molecule transition at smaller eb/ef- 
Indeed, we even find a transition if we just consider 
the bare molecule |A/o) and the "bare" polaron state 
|Po) — cJj^jiV), where particle- hole excitations of the 
Fermi sea are completely neglected. This "bare" bind- 
ing transition, in fact, corresponds to a spinodal in the 
mean-field theory of a two-component Fermi gas [isj in 
the limit of full polarization. To see this, we note that 
the equation for the bare molecule \Mq'^^): 



g ^ ^ Em - Ckt - eq-ki 

is equivalent to the linearized mean-field gap equation 
for the superfluid order parameter at momentum Q, with 
spin- up chemical potential fi^ — ef and spin-down chem- 
ical potential = E. At the binding transition E = 0, 
the normal unpaired phase becomes linearly unstable to 
forming a paired Fermi superfluid. For Q = 0, the tran- 
sition occurs at Eb/ef = r, as plotted in Fig. [2j Note 
that, for an infinitely massive impurity (r = 0), this im- 
plies that there is always a bound molecule when > 0, 
consistent with the exact solution ([5]). 

With the inclusion of particle-hole excitations, we find 
that the polaron-molecule transition survives provided 
we treat the polaron and molecule on an equal footing 
and consider wave functions \P) and \M). The effect 
of a single particle-hole pair is to shift the transition to 
higher eb /ef ior r > 0, but we see in Fig. [5] that the be- 
havior of the \P)-\M) transition line is still qualitatively 
similar to that of the bare transition line. The q approx- 
imation gives a transition line that is closer to the full 
\P)-\M) transition line than to the bare one, but it still 
underestimates eb/ef for a given r > and it suggests 
that the variational approach is least quantitatively ac- 
curate for the transition at r = 1. We find that the 
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FIG. 2: (Color online) Polaron-molecule transitions that have 
been determined using three different variational approaches 
of increasing accuracy for the polaron and molecule. To the 
left (right) of each transition line, the polaron (molecule) is 
the ground state. The "bare" transition line was determined 
using the bare variational wave functions |Po), \M^^-°^) and 
corresponds to a mean-field spinodal line (see text). The full 
calculation for the "dressed" variational wave functions |P), 
\M) yields the solid (black) line, while an approximate calcu- 
lation that averages the dependence of the energy on the hole 
momentum q gives the dashed (blue) line. 

variational wave function \P) is less accurate than \AI) 
at the transition and thus we expect the exact transition 
to lie at slightly higher Eb/ef for each r > 0. 

Having established the existence of a polaron-molecule 
transition in 2D, it is natural to ask whether there is a 
similar scenario in ID. The exact Bethe ansatz solution 
implies that there is no transition for r = 1 [l^. In- 
deed, we can prove that there are no polaron-molecule 
transitions for general r owing to the form of the sin- 
gularities in the variational equations in ID. Consider, 
first, the energy equation ([9]) for jA/g*^'), with the sum 
over k replaced by an integral: ^ X^k ~^ / {27t)<' ■ 
Q = IQI = fci?, we see that the integrand has a singular- 
ity at the energy of the bare polaron E = Em — £f = 
when k = Q. This singularity is integrablc for dimen- 
sions d ~ 2,S and so the interaction g is well-defined at 
this energy. However, for d = 1, we have a divergent 
integral and thus the energy of \Mq) can never equal the 
energy of \Po) for nonzero g or Eb, which means that the 
molecule never undergoes a transition to a polaron in this 
approximation. If we now consider the dressed polaron 
at finite momentum p: 
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then Eq. (|4]) becomes 
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where 



-E + Ep+q-ki + Ekt - Cqi- Hcrc, there is a 
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singularity in the q integral at the energy for |ii-jq 
when q = —kpp, which once again gives a divergent inte- 
gral in ID. Thus, we have E — > — oo, i.e. Sb ^ oo, as the 
energy of |P) approaches that of |Mo) from below. This 
singular structure is also present in higher-order varia- 
tional wave functions in ID, so that the higher-order 
molecule \M) always has a lower energy than \P) and 
so on. The absence of a transition in ID is perhaps be- 
cause there is no clear distinction between the statistics 
of a polaron and molecule in ID, and so we can expect a 
crossover from the weak to strong-coupling regimes. By 
contrast, in higher dimensions, a dilute gas of impuri- 
ties will exhibit different behavior (Fermi liquid or Bose 
superfluid) depending on whether they form polarons or 
molecules. 

Discussion A remaining question is whether or not 
the dressed molecule can have a finite momentum Q in 
the ground state for r > 0. It has been shown that this 
can give rise to a spatially-modulated superfluid for a low 
density of fermionic impurities [l7j . In the bare case, once 
£b/£f < ?'(r+l), we find that the molecule |Afo) acquires 
a finite momentum Q = y^kpa^/r — 1, where we have 
defined the lengthscale a = l/VSm^es using the reduced 
mass TTLr = i'T^i^ + m^^)^^- Thus, we find that the 
binding transition is actually at eb/sf = + with 
Q = kp at this point, in contrast to the 3D case where we 
always obtain Q = in the regime r < 1. The physical 
interpretation in 2D is that the bare molecule binds when 
Eb equals the center-of-mass kinetic energy of a molecule 
at Q = kp. However, further work is needed to determine 
whether or not this also occurs in the case of the dressed 
molecule and polaron. 

The molecule momentum Q at the polaron-molecule 
transition can affect the sharpness of the transition. For 
Q = 0, we require the creation of an extra particle-hole 
pair at the Fermi surface when the molecule decays into 
a polaron (or vice versa) in order to conserve energy 
and momentum. Thus, we expect to have a "first-order" 
transition where we can have a metastable molecule (or 
polaron) beyond the transition 0,0, HO]. However, for 
Q = kp ai the transition, we do not require the creation 
of a particle-hole pair and thus the transition is continu- 
ous, with an excited molecule (polaron) decaying quickly 
into a polaron (molecule) beyond the transition jl2l |. 

To connect with real 2D systems, we would ultimately 
like to know how the behavior of a single impurity ex- 
tends to that of a finite density of spin-down particles. 
In particular, is the polaron-molecule transition thermo- 
dynamically stable? Clearly, this will depend on the 
statistics of the impurity. For a two-component Fermi 
gas, a mean-field analysis 18| of Eq. ^ shows that the 
bare polaron-molecule transition is preempted by a first- 
order phase transition from an unpaired normal phase to 
a paired superfluid at Eb/sf = vT+T — 1. However, 



mean-field theory is unlikely to give reliable results in 
2D, so further work is required to better estimate the en- 
ergy of the superfluid and/or determine the interactions 
between dressed impurities. In principle, atomic gas ex- 
periments can probe the polaron-molecule transition by 
measuring |q;o P for the polaron using RF spectroscopy [3| 
or by determining the effective mass of the impurity using 
low- lying compression modes of the gas 0] ■ 

To conclude, we have shown that the inclusion of a sin- 
gle particle-hole excitation in the variational wave func- 
tion for the molecule is sufficient for reproducing the ex- 
act energy of an infinitely-massive impurity in the strong- 
coupling limit. Moreover, by using variational wave func- 
tions that capture both the weak and strong-coupling 
regimes, we have demonstrated that polaron-molecule 
transitions do exist in 2D but not in ID. 
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